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Abstract 

In order to test if the universal aspects of Barkhausen noise in magnetic 

materials can be predicted from recent variants of the non-equilibrium zero 

temperature Random Field Ising Model (RFIM), we perform a quantitative 

study of the universal scaling function derived from the Barkhausen pulse 

shape in simulations and experiment. Through data collapses and scaling 

relations we determine the critical exponents r and 1/avz in both simulation 

and experiment. Although we find agreement in the critical exponents, we find 

differences between theoretical and experimental pulse shape scaling functions 

as well as between different experiments. 
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I. INTRODUCTION 



Real materials, have dirt, or disorder, which often leads to slow (glassy) dynamics due 
to complex free energy landscapes with diverging energy barriers On long length scales 
and practical time scales, thermal effects often become unimportant. Indeed, when such a 
system is driven by an external field it jumps from one local metastable free energy minimum 
to the next, and the state of the system depends on its history. 

In magnetic materials jumps from one local minimum to the next involve a collective 
process whereby clusters of magnetic domains change the direction of their magnetization, 
in an avalanche. These avalanches can be triggered by a slowly but continuously increasing 
homogeneous external magnetic field H (taken from — oo to +00). These avalanches produce 
so-called Barkhausen noise, which can be observed experimentally as a voltage signal induced 
in a pickup coil wound around the magnet. Experiments show that these avalanches come in 
all sizes; their sizes are typically distributed according to a power law over several decades. 
Other systems also exhibit a broad range of avalanche sizes and durations following power 
laws: superconducting vortex line avalanches , resistance avalanches in a superconducting 
film ||, capillary condensation of helium in Nuclepore acoustic emissions in athermal 
martensites ||, earthquakes ||, and many others 0. 

We examine two variants of the non-equilibrium zero temperature random field Ising 
model (RFIM), which both predict power law distributions of avalanche sizes and universal 
non-equilibrium collective behavior. The RFIM is a model for a conventional magnet in 
which magnetic domains are modeled by "spins" on a lattice that can only point up or 
down. The first variant basically exhibits single front (domain wall) propagation dynamics 
in which spins at the edge of an existing front flip when it is (locally) energetically favorable 
to do so. Spins that are not adjacent to this front are very unlikely to flip on their own 
due to the presence of an infinite-range demagnetizing field, which is present in addition to 
nearest neighbor ferromagnetic interactions. We call this the front propagation model ||. 
In the second variant, called the nucleation model, spins flip anywhere in the system when 
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it is (locally) energetically favorable to do so. In this case there are many interacting fronts, 
unlike in the front propagation model. Only nearest neighbor ferromagnet interactions are 
included in the nucleation model ||. 

In this paper we perform a detailed quantitative comparison of the universal avalanche 
pulse shapes and exponents obtained from: front propagation dynamics, domain nucleation 
dynamics, and experiment. Our analysis constitutes a test of whether the non-equilibrium 
zero-temperature RFIM (either variant) is in the same universality class as experimental 
systems exhibiting Barkhausen noise. 

II. THE MODEL 

The RFIM consists of a (hypercubic) lattice of iV spins (sj = ±1), which may point 
up (sj = +1) or down (sj = —1). Spins are coupled to nearest neighbors (through a 
ferromagnetic exchange interaction J), and to an external field H(t) which is increased 
adiabatically slowly. To model dirt in the material, we assign a random field, hi, to each 
spin, chosen from a distribution P(hi) = exp(—h^/2R 2 )/^/2irR, where R, the disorder, 
determines the width of the Gaussian probability distribution and therefore gives a measure 
of the amount of quenched disorder for the system. 

The Hamiltonian for the system at a time t is given by: H — —JsiSj — J2i{H(t) + hi — 

<ij> 

J in fM)si, where Jj n / is the strength of the infinite range demagnetizing field (Jj n / = for the 
nucleation model), M = s % is the magnetization of the system, and < ij > stands for 

i 

nearest neighbor pairs of spins. Initially, H(—oo) = —oo and all the spins are pointing down. 
Each spin is always aligned with its local effective field h* = JJ2<ij> Sj + H(t) + hi — J in fM. 
The external field Hit) is adiabatically slowly increased from — oo until the local field, hi , 
of any spin Si changes sign, causing the spin to flip [p|,|T0f. It takes some microscopic time 
At for a spin to flip (At = 1 for our simulation). The spin flip changes the local field of the 
nearest neighbors and may cause them to flip as well, etc. This avalanche process continues 
until no more spin flips are triggered. Each step of the avalanche, that is, each At in which 
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a set of spins simultaneously flip, is called a shell. The number of spins that flip in a shell 
is directly proportional to the voltage V(t) during the interval At that an experimentalist 
would measure in a pick-up coil wound around the sample. In our simulations we therefore 
denote the number of spins flipped in a shell at a time t by V(t). The first shell of an 
avalanche (one spin flip) is triggered by the external field H(t), while each subsequent shell 
within the avalanche is triggered only by the previous shell, since H(t) is kept constant 
while the avalanche is propagating. H(t) is only increased when the current avalanche has 
stopped, and is increased only until the next avalanche is triggered (i.e. — > 0). The 
number of shells in an avalanche times At defines the pulse duration, T, or the time it took 
for the entire avalanche to flip. In this paper we will be interested in looking at V(t,T) for 
< t < T, that is, the voltage as a function of time for an avalanche of a given duration T. 

The front propagation model exhibits self-organized criticality (SOC) fll"|-|l~3|| . This 
means that as H is increased the model always operates at the critical depinning point, 
and no parameters need to be tuned to exhibit critical scaling behavior (except ^ — > 0). 
The nucleation model, on the other hand, is a plain old critical system with a continuous 
second order phase transition with disorder as the tuning parameter. The continuous non- 
equilibrium phase transition can be understood as follows. For zero disorder, the random 
fields of each spin will be the same, so that when one spin is flipped the entire lattice of 
spins will flip. This results in a rectangular hysteresis curve with a macroscopic jump in the 
magnetization when the external field overcomes the interaction with the neighbors. On the 
other hand, if the disorder is infinite, each spin will have a very different random field, so as 
the external field is raised each spin will essentially flip independently, triggering no other 
spins to flip. This will result in a smooth "hysteresis" curve with an approximately constant 
slope (M(H) ~ H over a wide range of H). In between these two phases of behavior there is 
a continuous phase transition at some critical value R = R c . At the transition each branch 
of the hysteresis loop has a single point with infinite slope (^ \±h c — > oo) at a critical field 
H = ±H C . R c and H c are nonuniversal. In three dimensions they are: R c = 2.16 and 
H c = 1.43 (in units of J) [14]]. The transition is characterized by a number of universal 
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critical exponents [15| , and scaling functions. In this paper we specifically focus on only two 
exponents: 1/avz and r, and one scaling function. The significance of these exponents is 
given later. The scaling function we examine can be obtained through experiment as well 



as simulation. Details of the simulation algorithm are given elsewhere [|16[ 



III. THE EXPERIMENTS 

In addition to examining results obtained from simulation we study results from three 
different experiments. We performed one of these experiments, and the results for the other 
two experiments were obtained from already published results [|l7l ,|lS|1. 

Our own experiment was performed on an (unstressed) amorphous alloy, Fe2\Co§/s t B\*,. 
The data we present from Durin et al |18| is from an experiment on an amorphous alloy with 



a different composition, Fe 64: Co2iBi 5 , under a tensile stress. Studies of the effect of tensile 
stress on samples of this type indicate that for low tensile stress, the domain structure is 



a complicated pattern of maze domains, dominated by quenched-in stresses pJJ. On the 



other hand, when such materials undergo tensile stress, the uniaxial anisotropy gives way 



to a simpler domain structure with a few parallel domains in the direction of the stress |19 
Related to this change in domain structure is a change in the dominant interaction in the 
material. In amorphous alloys under stress, surface tension effects are thought to be more 
important than dipolar interactions, while dipolar interactions dominate for polycrystals 
and materials with small grains [pOfl . 

Durin et al's sample was under stress so as to enhance stress-induced anisotropy so much 
that the long range dipolar interactions can be neglected, placing their experiment into the 
universality class of the front propagation model [p0| , p2| , p3|1 . Spasojevic et al's sample was 
a quasi-two-dimensional metal glass, more precisely a commercial VITROVAC 6025 X [|17| . 



Based on the scaling exponents obtained from Spasojevic et al's experiment (given later), 
their experiment does not seem to fall into any universality class discussed in this paper. 
Our experiment seems to be in a crossover regime between two universality classes; details 
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are given later in this paper. Further details about Durin et al's, and Spasojevic et al's 
experiment can be obtained elsewhere [pj-PDf . 

IV. EXTRACTING PULSE SHAPES 

Nucleation Model: We simulate four realizations of a 1200 3 system near R = R c (R = 2.2) 
and record the time series V(t,T) of avalanches from an H window near H c (1.42 < H < 
1.43). We average avalanches of a fixed pulse duration T (within the interval [T, 1.05T]), for 
various values of T. In each case we average over 1000 to 2000 avalanches to ensure strong 
fluctuations have been averaged out. We check finite size effects by performing simulations 
at 800 3 and 1000 3 and verifying that identical avalanche shapes (within small fluctuations) 
are obtained for all system sizes. 

Front Propagation Model: We perform 100 realizations of a 400 3 system and record 
avalanches from an H window within the slated part of the hysteresis loop (1.25 < H < 1.88). 
Even though the front propagation model exhibits SOC due to the infinite range demagne- 
tizing interaction, in order to avoid effects due to initial nucleation of the front (beginning 
of the hysteresis loop), or when the front encounters the boundaries of the simulation (the 
end of the hysteresis loop) we must choose avalanches near the middle of the hysteresis loop. 
We obtain avalanche shapes V(t, T) in a manner identical to how they were obtained for the 
nucleation model. We check finite size effects by performing simulations of 200 3 and 300 3 
size systems, and find a consistent avalanche shape for all three system sizes. 

Experiment: Measurements were performed on a 21 cm x 1 cm x 30 /im ribbon of 
Fe2iCo64-Bi5 alloy, a soft amorphous ferromagnet obtained from Gianfranco Durin. The 
domain walls run parallel to the long axis of the material, with about 50 domains across the 
width. A solenoid, driven with a triangle wave, applies a magnetic field along the long axis 
of the sample. Since domain wall motion dominates over other means of magnetization in 
the linear region of the loop, data were collected in only a selected range of applied fields 
near the center of the loop. The Barkhausen noise was measured by a small pick-up coil 
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wound around the center of the sample. This voltage signal was amplified, anti-alias filtered 
and digitized, with care taken to avoid pick-up from ambient fields. Barkhausen noise was 
collected for both increasing and decreasing fields for 80 cycles of the applied field through a 
saturation hysteresis loop. The driving frequency was 0.01 Hz; this corresponds to c = 0.09, 
where c is a dimensionless parameter proportional to the applied field rate and is defined in 
the Alessandro Beatrice Bertotti Montorsi model (ABBM model) for the Barkhausen effect 



24j[ . In this way, our measurements should be well inside the c < 1 regime identified in the 



ABBM model, in which we can expect to find separable avalanches rather than continuous 
domain wall motion. 

Due to background noise in real experimental data there is no definitive way to determine 
when an avalanche begins and when one ends; we set a sensitivity threshold which is high 
enough to cut out background noise and low enough to capture enough of the avalanche so as 
not to affect the shape. We check the validity of our threshold by perturbing the threshold 
by a small amount and noticing that there is no change in the avalanche shape. 



V. CRITICAL EXPONENTS AND DATA COLLAPSES 

In the front propagation model, as well as in the nucleation model near (R C ,H C ), the 
voltage V(T,t) scales as ||: 

V(T,t) = T 1 /™ z - 1 f shape (t/T) (1) 

By collapsing average avalanche shapes of various durations T we determine the universal 
scaling function, f s h ape {t/T), and the critical exponent 1/avz. The exponent \jovz relates 
the avalanche size, S, to the avalanche pulse duration, T, at criticality by S ~ T 1 ^" 2 . 
We find that the critical exponent 1/avz for the front propagation and nucleation model, 
obtained from simulation, is in close agreement with previous theoretical predictions and 
with experimental values [P5fl . The collapses are shown in Figs. 1 and 2. 

We also determine the avalanche size distribution that scales as D(S) ~ S~ T at criticality. 
The avalanche size distributions for our simulation and our experiment are given in Fig. 3; 
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the values obtained for r, for the front propagation model and nucleation model, are in 
close agreement with previously quoted values P,|14||. Experimentally, the scaling exponent 
t weakly depends on c in some materials ||19|| . From this dependence we find that there may 
be a difference of 0.02 between the value of r we find in our experiment and the value of 
t at c = (zero frequency). However, this difference is within the error bars we give for 
r. Table I summarizes results for r and \ javz for experiment, simulation, and mean field 
theory. We see that the critical exponents for the front propagation model agree (within 
error bars) with exponents from Durin et al's experiment, as expected. 

While Durin et al's experiment is believed to fall into the front propagation universality 
class, as discussed above, our experiment is neither in the front propagation nor in the 
mean field universality class. A priori we would assume that our experiment would be in 
the mean field universality class since our sample was unstressed, and previous experiments 
have indicated that unstressed samples will exhibit mean field behavior [2(|. The critical 
exponents found from our experiment indicate that it maybe in a crossover regime between 
the mean field and the front propagation universality classes. The exponent r has a value 
of 1.46 ± 0.05 (see also Table I), which is between the value of 1.28 for the front propagation 
model (sample with stress) and the mean field value of 1.5 (sample without stress) P,|2"0" 



In addition, the exponent a = 1.74 ± 0.06, determined from the avalanche pulse duration 
distribution given by D(T) ~ T~ a at criticality, (see Inset of Fig. 2), is between the value 
of 1.5 for the front propagation model (with stress), and the mean field value of 2 (without 
stress) |18|]. Residual stress on our sample may have resulted in these anomalous exponents. 

Although our experiment may be in a crossover regime, we were able to obtain a good 
collapse of the avalanche pulse shapes (see Fig. 2). In order to reaffirm the validity of the 
exponent obtained from the collapse, we independently checked the value of 1/auz from 
the power spectra (P(w) ~ w~ l l (TUZ at criticality), and found l/avz = 1.73 ± 0.08, which 
is consistent (within errors bars) with \javz = 1.70 ± 0.05 obtained from the experimental 
avalanche pulse collapse. 



S 



VI. FITTING TO ORTHONORMAL POLYNOMIALS 



The mean field shape of avalanches in our model V MFT (t, T) is an inverted parabola 
| ZT]|. In order to study corrections to the mean field shape, in simulation and experiment, 



we derive a set of orthonormal polynomials fi(t) with /*(— L) = fi(L) = 0, where L = T/2 is 
half the duration of the avalanche. The negative of the first polynomial, fo(t), is proportional 
to the mean field result. The first four polynomials of the set are: 



Mt) = J^P-L>) (2) 



/2(t) = V6i^ (7t4 ~ 8LV + L4) (4) 



/3(t) = V6^ (3t5 - 4L¥ + L4t) (5) 



h{t) = V 2^ (33t6 " 51L¥ + 19LH2 ~ (6) 
We fit a linear combination of the above polynomials to the average avalanche shape 
obtained from simulation and experiment. L is also left as a free parameter in the fits to 
the average shapes. The pulse duration is then precisely defined as T = 2L. The results 
of the fits are shown in Fig. 1 and Fig. 2. In Fig. 4 we give the coefficients for the fits, 
found in simulation and experiment. We also include the coefficients for the fit to an average 



avalanche shape determined by Durin et al's experiment [18], whose avalanche shape is given 



in Fig. 5. The coefficients are determined from the total fit function: 

F(t) = aMt) + ai/i(t) + a 2 / 2 (t) + a 3 f 3 (t) + a 4 f 4 (t) (7) 

a , ci2, and a 4 are the symmetric coefficients of the fit, while a 1; a 3 are the antisymmetric 
coefficients of the fit (i.e. multiplying polynomials that are not symmetric under time reversal 
t — > —t). We are particularly interested in the asymmetry of the avalanche shapes. 

From inspection of the avalanche shapes we see that the experimental avalanche shapes 
are strongly asymmetric under time reversal, while the avalanche shapes determined from 
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the simulation of the two models are both very close to symmetric. Quantitatively, the 
coefficients for the avalanche shapes of the nucleation model are very similar to those of 
the front propagation model and both are different from the experimental avalanche shapes' 
coefficients. While avalanche shapes in both models are slightly asymmetric to the left 
(i.e. V(t,T) increases more slowly than it decreases), the experimental avalanche shapes 
are strongly asymmetric to the right direction (i.e. V(t,T) increases quickly and decreases 
slowly). The origin of this difference between theory and experiment is not yet understood. 
It comes as a surprise, since the critical exponents obtained in Durin et al's experiment are 
in close agreement with exponents obtained from the front propagation model (see Table I) 
and the avalanche shapes are expected to be just as universal as the critical exponents. On 
the other hand, we see in Fig. 5 that different experiments give different pulse shape scaling 
functions. This difference may be a result of the fact that the experiments are not in the 
same universality class. Nevertheless, all the experimental pulse shapes do have the same 
sign of asymmetry. 

We did one check on whether the pulse shape function in our model is indeed universal, 
an assumption often taken for granted: We changed the lattice of our simulation from a 
simple cubic to a BCC (body centered cubic) lattice and found an identical pulse shape 
scaling function (within small statistical error), as expected when universality holds. 

VII. DISCUSSION AND CONCLUDING REMARKS 

In performing these analyses we have stumbled upon an interesting observation: after 
appropriate rescaling of the y-axis, the universal pulse shape function for the front prop- 
agation and nucleation model appear to look the same (see Inset for Fig. 1) even though 
we know J7| that the models do not belong to the same universality class. However, upon 
more precise quantitative analysis we see that they are in fact different, their a\ coefficients 
are different by several a. Furthermore, we tried to collapse two pulse shape functions with 
T ~ 52 from the two different models and found that they could not be collapsed precisely 
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despite, naturally, scaling by different critical exponents appropriate for the two models. The 
front propagation pulse scaling function is more asymmetric, as supported by the somewhat 
larger a\ fit coefficient compared to the nucleation model result. 

By examining the pulse shape scaling functions as a sharper test than merely critical 
exponents for the universality class of the non-equilibrium zero temperature RFIM, we 
raised many questions. What accounts for the difference between theory and experiment, 
and between different experiments? Is the theory incomplete or inaccurate at this level of 
description? Experimentally, we do not yet know what material features are required to 
produce universality of the pulse shape function. Differences between experimental results 
make further experimental tests desirable (see Fig. 5). 

Although the three experiments we examined do not fall into the same universality class 
and do not have very similar pulse shapes, they share one universal aspect not shared by 
the results from simulation: all the experimental averaged pulse shapes are asymmetric in 
the leftward direction. This suggests that there may be a phenomena that exists in real 
materials that has not been accounted for in theory. 
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FIGURES 
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FIG. 1. Pulse shape collapses of V(t,T) obtained from simulation. Three pulse shapes were 
collapsed for each model; these pulse shapes represent averaged avalanches of pulse durations 
T = 52, 73, and 106 within 5%. From the collapse of front propagation model avalanche pulse 
shapes, we find \ jovz = 1.72 ±0.03. From the collapse of nucleation model avalanche pulse shapes, 
we find 1/avz = 1.75 ± 0.03. The bold line going through the collapses is the non-linear curve 
fit obtained from the set of orthonormal polynomials presented in this paper (Eqns.(2)-(6)). Note 
that the non-linear curve fit is shown for only one of the collapsed averaged avalanches in each 
case. Inset: By rescaling the height of the nucleation model collapse by 20% we obtain a collapse 
of pulse shapes from the two different models suggesting that their pulse shapes are very similar, 
but quantitatively not the same as described in the text. 



12 




FIG. 2. Collapse of averaged experimental avalanche pulse shapes is shown, yielding 
1/avz = 1.70 ± 0.05. The four curves represent averaged avalanches of pulse duration: T = 88t a , 
110i o , 132£ , and 165i D within 10%, where t Q = represents the time between each measure- 

ment of the Barkhausen noise train. Each of the four curves is an average of between 1152 to 
1561 avalanches. The smooth bold curve is a fit of the averaged avalanche of duration T = 132i Q 
using the orthonormal polynomials given in Eqns.(2)-(6). Inset: Distribution of avalanche pulse 
durations, D(T), obtained from our experiment. Avalanche pulse shapes were extracted from the 
region indicated by the arrows, and this region is well within the scaling regime. In this scaling 
regime D(T) scales as D{T) ~ J 1 -" where a = 1.74 ± 0.06. This value of a is between the mean 
field value of a = 2.0 (sample without stress) and the value of a = 1.5 for front propagation 
dynamics (sample with stress) pPf , indicating that our experiment may be in a crossover regime. 
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FIG. 3. Avalanche size distribution for front propagation model, nucleation model, and our 
experiment. The exponent r given in the legend is the critical exponent corresponding to the 
scaling of the avalanche size distribution (D(S) ~ S~ T ). 
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FIG. 4. Fitting coefficients to the avalanche shapes determined for the two models, our 
experiment, and Durin et al's experiment [18|. We find that the coefficients are very similar 



for the two models. While the a± coefficients determined from the experiments are significantly 
different from the two models, the difference is not only in the sign of asymmetry but also in the 
magnitude of asymmetry. Each fitting coefficient, except for Durin et al's, was determined from 
three realizations of the universal scaling function in each case. The coefficients, for the two models 
and our experiment, plotted above, represent median values, while the error bars are determined 
from the higher and lower values. Durin et al's avalanche shape, presented in Fig. 5, was used to 
calculate the coefficients presented above; no error bars are provided in this case since only one 
realization was available. 
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FIG. 5. Comparison of our experimental pulse shapes with experimental pulse shapes obtained 
by two other groups. Durin et al's sample was under stress with large stress induced anisotropies, 
putting their experiment into the universality class of front propag ation For Durin's 

experiment 1/avz = 1.77 [18| and r = 1.27 |20j,22|. For Spasojevic's experiment 1/avz = 1.58 and 
r = 1.77 0. 
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TABLES 





Nucleation 
Model 


Front Propagation 
Model 


Our Experiment 
Durin et al's Expt.t 
Spasojevic et al's Expt.* 


Mean Field 1 16] 


T 


1.60 ±0.04 
1.60 ±0.06 \14\ 


1.33 ±0.08 
1.28 @ 


1.46 ±0.05 
1.27±0.03 t |2C,22] 
1.77* 


1.5 


1/avz 


1.75 ±0.03 
1.75 ±0.07 [14] 


1.72 ±0.03 
1.72 


1.70 ±0.05 
1.77±0.12t |2C,18] 
1.58* |TJ 


2 



TABLE I. In the above table we present the critical exponents r and 1 jovz in d = 3 dimensions 
for the nucleation model, front propagation model, and all three experiments discussed in this 
paper. We also include the mean field values of these exponents. 
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